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ABSTRACT 

'Deep surveys at 850/im from Mauna Kea using the SCUBA camera on the JCMT appear to have discovered a 
5_i -substantial population of ultraluminous infrared galaxies (ULIGs: Ljr > lO^^L0)0at high redshift {z ~ 1—4). 
CL|The cumulative space density of these sources (~ 10^ deg^^ with 5*850 > 1 mJy) is sufficient to account for 
-nearly all of the extragalactic background light at submillimeter wavelengths. Current estimates of the 
^ redshift distribution suggest a peak in the comoving space density of SCUBA sources at ^; = 1 — 3, similar 
to what is observed for QSOs and radio galaxies. The luminosity density in the far-infrared/submillimeter 
.exceeds that in the UV by factors of 3-10 over this redshift range, implying that as much as 80-90 % of the 
^—1 "activity" in galaxies at z 4 is hidden by dust. The SCUBA sources plausibly represent the primary epoch 
in the formation of spheroids and massive black holes triggered by major mergers of gas-rich disks. 
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INTRODUCTION 



The Submillimeter Common User Bolometer Array (SCUBA) camera (Holland et al. 1999) recently installed 
on the James Clerk Maxwell Telescope (JCMT) has provided a new window for ground-based studies of the 
|high-z Universe. This brief review summarizes results from the large campaign of deep submillimeter surveys 
ii.carried out during SCUBA's first year of operation on Mauna Kea (~ Fall '97-Summer '98). Evidence is 
O presented that the SCUBA detections must be predominantly ULIGs at high redshift, and that the cumulative 
^ -space density of these sources is sufficient to account for nearly 100 % of the extragalactic background light 
^ (EBL) at far-infrared/submillimeter wavelengths. Progress as well as current problems in identifying optical 
counterparts to the SCUBA sources are then discussed, followed by a review of model calculations comparing 
.5^ ]the luminosity density in the submm/far-IR with that in the UV over the redshift range, z = — 5. Finally, 
^ -we discuss evidence which suggests that the SCUBA sources, like local ULIGs, are powered by intense 
^ .circumnuclear starbursts and powerful AGN both of which are fueled by major mergers of gas-rich disks. 

SCUBA DEEP SURVEYS 



As this Workshop was being held, several observing teams were announcing results from the first deep 
blank- field submillimeter surveys at 450;um/850/im using SCUBA. Papers by Hughes et al. (1998) for the 
UK-consortium and by Barger et al. (1998) for the Hawaii-Japan group were announced in press releases 
during the week of the workshop, and results from the Canada- France consortium were about to be submitted 
(Eales et al. 1998). These "ultradeep" and "deep" blank-field surveys follow the pioneering surveys of Small et 
al. (1997) who were the first to infer a substantial population of luminous submillimeter galaxies from their 
SCUBA detections at 850 /im of background sources amplified by weak lensing from foreground clusters. 
With the exception of one source detected at 450 //m in the cluster fields, all of the SCUBA deep field 



^Invited review for the workshop "Space Infrared Telescopes and Related Science", at the 32nd COSPAR Meeting, Nagoya, Japan, 
1998, eds. T. Matsumoto and T. de Graauw. To be published in Advances in Space Research (Oxford: Elsevier), 1999. 

^This paper draws on new material from the publication by Barger et al. 1999, and the website http:/www. ifa.hawaii.edu/~cowie. 
Due to previous commitments, Len Cowie and Amy Barger were unable to present this review. 

^Lir = L(8 — 1000) /im in the object rest-frame. Unless stated otherwise. Ho ~ 75kms~^Mpc~^, qo — Q . 



detections have been at 850/im. The HDF (Hughes et al. 1998), and Lockman Hole + SSA13 fields (Barger 
et al. 1998a, b) are currently the three fields with the deepest SCUBA observations, each with ~ 50 Hrs total 
integration time reaching a noise level, cisso ~ 0.6-0.8 mJy. Our results for the Lockman Hole are discussed 
below. 



The Lockman Hole 



The 850 /xm data on the Lockman Hole field (LH_NW: Barger et al. 1998, 1999) are shown in Figure 1. The 
two SCUBA sources, LH_NW1 and LH_NW2, detected at 850 ^um (> 3 a), have 850 /xm fluxes of 5.1 mJy, and 
2.7 mJy, respectively, with upper limits at 450 //m of ^50mJy (Sex). 




Figure 1. SCUBA 850 /im detections (smah thick circles: Barger et al. 1998, 1999), and ISOCAM 
7/xm detections (small thin circles: Taniguchi et al. 1997) in the Lockman Hole northwest (LH_NW) 
Deep Field (J2000: RA = 10'*33"'55.5^ Dec = +57°46'18") superimposed on a K' image obtained 
with the QUick Infrared Camera (QUIRC) on the University of Hawaii 88-inch telescope. The 
field-of-view of the ISOCAM detector array and the SCUBA array are indicated by a long dashed 
line and large solid circle respectively. On the right are two zoomed images of the region outlined 
by the 45" x 45" box which is centered on the strongest SCUBA source The zoomed K' image was 
obtained with the Near InfraRed Camera (NIRC) on the Keck 10-m telescope, and the zoomed 
B-band image was obtained with the University of Hawaii 2.2-m telescope. 

Deep observations of the LH_NW field have also been obtained in the mid- and far-infrared with ISO, at B 
and K' from Mauna Kea, and at 20 cm with the VLA. As shown in Figure 1, neither SCUBA source has an 
ISOCAM 7/im counterpart (< 35 /uJy; 5(t). LH_NW1 appears to be centered on a faint K' source (i^^B ~ 
21.8) with disturbed morphology, which is barely detected in the current B-band image (-Bab = 23.5). 
LHJNrW2 is "blank" in our K' and B images implying that any counterpart has > 22.5 and Bxb > 24.5. 
The ISOPHOT sources detected in a wider area (44' x 44') survey centered on LH_NW have relatively large 
position uncertainties (typically 20"-30"), however there are no 95 /xm or 175 /xm detections within ~2' of 



either SCUBA source (the 5 a upper hmits are Sg^ < 40mJy and 5*175 < 75mJy: Kawara et al. 1998). The 
VLA results are still preliminary, but any radio continuum counterparts must have 5200111 ^ 500 /xJy (5 a) 
(M. Yun, private communication). 

ULIGs at High Redshift 

Prom the strength of the 850 jim detections and the faintness of the K' counterparts alone, it is relatively 
straightforward to show that the SCUBA sources are most likely to be ULIGs at high redshift (i.e. z >1). 
The "submillimctcr excess", (= uS,,{ShQ) /vS^{2.2)), for both LH_NW1 and LH_NW2 is larger than 1 (2.4 
and > 3 respectively), which is impossible to produce from normal optically selected galaxies at any redshift, 
or even by the most extreme infrared selected galaxies at low redshift, but is almost exactly what would be 
expected for an ULIG at high redshift. Figure 2 shows that the expected flux for the prototypical ULIG 
Arp 220 at z > 1 is on the order of a few mJy at 850 /xm. Also, the combination of a large negative K- 
correction in the submillimeter plus a relatively flat or positive K-correction in the near-IR naturally leads to 
values vSy{S^Q)/uSu{2.2) > 1 for all ULIGs at z ^1.5. The observed faintness of the high-z submillimeter 
sources in current B-band images and the non-detections at 7/Ltm in the deep ISOCAM images are consistent 
with the large U-B colors and the pronounced minimum at ~ 3-6 /um respectively, in the rest-frame SEDs of 
ULIGs like Arp 220. The non-detections at 95 /um and 175 ^m are consistent with the relatively high noise 
levels in the ISOPHOT images coupled with the steep positive K-correction encountered over the rest-frame 
wavelength range ~ 20-80 /Ltm for all ULIGs. 
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Pigure 2. Observed radio-to- UV spectral energy distribution of the nearest ULIG, Arp 220 [z = 
0.018). Labeled tickmarks represent object rest-frame emission that will be shifted into the 850 fim 
and 2.2 jim observed frame for redshifts, z = — 5. The insert shows the corresponding observed- 
frame 850 /xm flux and vSui850)/i'Su{2.2) ratio for Arp 220 at redshifts, z = - 5. 



SOURCE COUNTS AT 850//m 



Table 1 summarizes the 850 fim cmnulative source counts from all of the SCUBA deep field surveys published 
to date. Taking into account the different sensitivity limits and the reported formal uncertainties, there 
appears to be good agreement on both the shape of the logA^'-logS plot over the observed range of 2-10 mJy, 
as well as on the extrapolated cumulative counts down to 1 mJy obtained from P{D) analyses of the low-level 
fluctuations in the ultradeep images. In retrospect, it is remarkable how little time it has taken to reach 
consensus on the surface density of faint submillimeter sources per SCUBA field. This success is first and 
foremost due to the excellent performance of SCUBA, but was also helped by a phenomenally well timed six 
consecutive months of near record low mean submillimeter opacity on Mauna Kea, and the commitment of 
relatively large blocks of observing time to several teams of observers. 



Table 1. SCUBA Source Counts at 850 microns' 
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2 clusters (A370, CL2244-02) 

background sources in 3 star fields 

LH_NW; SSA13 

HDF 

CFRS 

7 clusters 

LH_NW; SSA13; SSA17; SSA22 



modified and updated from a compilation by Trentham et al. (1998) 



850 2500±1400 4 

1100±600 8 

800l^J°° 3 

7000±3000 1 

1800±600 2.8 

1500±700 4 

7900±3000 1 

10,000±2000 1 



Small et al. (1997) 
Holland et al. (1998) 
Barger et al. (1998) 
Hughes et al. (1998) 
Eales et al. (1998) 
Blain et al. (1999) 

Barger et al. (1999) 
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Figure 3. Comparison of the 850 //m source counts (solid squares: from Barger et al. 1999) with 
semi-analytic model counts (see text). Model 1 is equivalent to the "no evolution" model of Blain 
&: Longair (1996). Models 2 and 3 are equivalent to Models A and E respectively from Guiderdoni 
et al. (1998). 

Figure 3 shows that the cumulative 850 /xm counts from 2 mJy to 10 mJy can be approximated by a single 
power law of the form N{> 5) = 1 x 10^ S~'^ deg~^. The extrapolated value at 1 mJy is consistent with the 
P{D) results reported in Table 1. This large surface density implies strong evolution in the ULIG population. 



Figure 3 compares the observed SCUBA counts with predictions from semi-analytic models using three rather 
extreme distributions of ULIGs. Model 1 is based on the local IRAS 60 ^um luminosity function of galaxies 
(i.e. ~ 0.001 ULIGs deg-^ at z < 0.08: Kim & Sanders 1998; see also Soifer et al. 1987; Saunders et al. 1990) 
assuming no evolution, which underestimates the observed sources by nearly 3 orders of magnitude. Model 
2 includes no ULIGs, but instead attempts to account for the fraction of the optical/UV emission absorbed 
and reradiated by dust in sources observed in optical/near-UV deep fields. Model 2 still substantially 
underpredicts the 850 //m source counts by a factor of ~ 30. Model 3 includes a strongly evolving population 
of ULIGs, constrained only by recent measurements of the submillimeter background. Model 3 is a much 
better match to the observed SCUBA counts, although the particular redshift distribution assumed (Model 
E of Guiderdoni et al. 1998) was arbitrarily chosen to push most of the ULIGs to z > 3, which as described 
below is at odds with current redshift estimates. 

CONTRIBUTION OF SCUBA SOURCES TO THE SUBMM/FAR-IR BACKGROUND 

Figure 4 compares the contribution of the 850 fim sources with the recent model of the EBL determined 
from COBE data at far-infrared/submillimeter wavelengths (Fixsen et al. 1998; see also Hauser et al. 1998; 
Puget et al. 1996, 1998). Approximately 25% of the 850 /^m background resides in sources brighter than 
2mJy, and nearly all of the EBL at 850 nm can be accounted for by sources brighter than 1 mJy, assuming 
the extrapolation down to 1 mJy given by the straight line fit to the SCUBA data in Figure 3. 
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Figure 4. Comparison of the contribution of the 850 /xm sources brighter than 3 mJy (solid circle) 
and extrapolated contribution of sources brighter than 1 mJy (open circle) to the EBL compared 
with the Fixsen et al. (1998) analytic approximation (solid curve) to the EBL. The two dashed 
curves are for observed source temperatures of 50 K and 25 K where each is based on a A-weighted 
Planck function, (see [ittp: / /www.ifa.hawaii.edu/^cowie/| ) 

Figure 4 makes the important additional point that the average observed temperature of the 850 /im sources 
must be low (i.e. ^ 25 K) so that their emission will not exceed the EBL at wavelengths shortward of 
200 /im. If the SCUBA sources are indeed distant analogs of local ULIGs whose typical dust temperatures 
are ~ 50-80 K, then a low observed temperature would provide additional evidence that the SCUBA sources 
lie predominantly at high redshift (i.e. z > 1). 



SOURCE IDENTIFICATION 



Progress in identifying optical/near-IR counterparts of the SCUBA deep-field sources has been frustratingly 
slow, due in large part to the uncertainty in the 850 //m SCUBA positions, but also to the intrinsic faintness of 
many of the sources, even at K-band. A prime example is the current controversy surrounding the identified 
SCUBA counterparts in the HDF. Four 4 of the 5 850 ^m SCUBA sources in the HDF, the identifications by 
Hughes et al. (1998) hinge in large part on plausibility arguments using "template SEDs". Although these 
attempts are admirable, they apparently can be misleading, for example the identification of the brightest 
source, HDF850.1, with a faint counterpart at 2; = 3.4 is now known to be in error from a more accurate 
position obtained from millimeter interferometry. Our results for the Lockman Hole are not much more 
encouraging (see Figure 1): LH_NW1 seems to have an identified K-band counterpart, however no emission 
lines are seen in a deep Keck spectrum, and LH_NW2 is a blank field in the current K-band image. 

Progress has been somewhat better in the lensed cluster fields, where the 1-2 mags of amplification appears 
to reduce the number of blank fields (^ 20% at / < 25: Small et al. 1998), but uncertainty in the 850 //m 
positions is still a major problem, plus contamination by cluster members, unless properly accounted for, 
will artificially inflate the number counts and skew the distribution of redshifts. The wider-area shallower 
surveys (e.g. Eales et al. 1998; Lilly et al. 1999) preferentially select brighter nearby sources, which are more 
likely to have brighter optical/near-IR counterparts, but still suffer from position uncertainty at 850 /xm. 

The net result after combining all current studies is that only ~ 20-25% of the 850 //m sources currently have 
"secure" identifications, and even for these it has been necessary to rely heavily on data at other wavelengths 
to make a strong case that only the identified source has a high probability of being correct. 

TEMPLATE SEDs 

Attempts to sort out which of typically 2-3 equally probable sources is more likely to be the true SCUBA 
counterpart have made use of "template SEDs" to determine which source is most likely to match the known 
properties of ULIGs. This procedure works best when additional deep-field data are available at several 
wavelengths, usually deep radio continuum images from the VLA, and high resolution, deep near-IR and 
optical images from either HST or with adaptive optics from the ground. The most complete and deepest 
data sets are of the HDF, with shallower multiwavelength data generally available for other SCUBA fields. 
However, even for the HDF, measurement uncertainties coupled with a real dispersion in ULIG properties 
can still give ambiguous results. 

Figures 5 and 6 illustrate the range of rest- frame SEDs observed for a sample of low-redshift ULIGs, and the 
corresponding flux ratios that would be observed for these objects at high redshift. It is instructive to briefly 
review what has (or hasn't) been learned from trying to identify the correct SCUBA counterpart and/or to 
set meaningful constraints on redshift using observed flux ratios. 

Radio /Far-Infrared Correlation 

Perhaps the most promising means of identifying SCUBA positions is from high resolution radio continuum 
maps, as has been done quite successfully for lower redshift IRAS sources (e.g. Condon et al. 1990). The 
subarcsecond radio positions coupled with a relatively tight correlation between the radio and the far-infrared 
emission for infrared galaxies (e.g. Condon et al. 1991) could provide relatively unambiguous identifications 
of the current SCUBA deep field sources, provided that it is possible to reach a point source sensitivity of 
^ 10 /xJy at 3cm or ^ 30 /xJy at 20cm. Unfortunately, although this has been achieved for the HDF (Richards 
et al. 1998), the majority of the SCUBA sources in the HDF are not detected in the radio continuum at a 
level which is already at least 3-5 times fainter than even the maximum observed ratios shown on the plot 
in Figure 6. One possible explanation would be to push the majority of the SCUBA sources to very large 
redshifts, i.e. z > 3 (e.g. Carilli and Yun 1999), but as shown below, this is in conflict with other best 
estimates of the range of redshift for the majority of the SCUBA sources. The current situation is bizarre 
enough to have caused others to attempt interesting explanations such as a substantial underestimate of the 
pointing error in the JCMT observations of the HDF (Richards 1999). A more likely explanation is that the 
SCUBA sources are indeed fainter in the radio than ULIGs studied by IRAS, although why this should be is 
currently not clear. The bottom line is that ultradeep VLA observations are currently incapable of detecting 




Figure 5. Measured rest-frame spectral energy distributions (SEDs) from UV-to-radio wave- 
lengths for a representative sample of the most luminous sources detected by IRAS. 




Figure 6. Left panels: (top) flux ratios as a function of redshift for nearby ULIGs from the IRAS 
Bright Galaxy Sample (Soifer et al. 1987; Sanders et al. 1995); (bottom) Flux versus redshift for 
Arp220. Right panels: (top) B-R versus redshift (bottom) B-K versus redshift for three nearby 
ULIGs (dashed/dotted lines) from the IRAS Bright Galaxy Sample compared with the redshifted 
colors of a standard elliptical galaxy (Trentham et al. 1999). 



the majority of SCUBA deep field sources, and thus even deeper ultradeep VLA observations are needed if 
there is to be any hope of using this potentially powerful diagnostic to pin down the SCUBA positions. 

450/Ltm/850/um Spectral Index 

In all of the SCUBA deep-fields, only one source has been reliably detected (5cr) at 450 /xm, that being 
SMM02399-0136 in the field of the z = 0.37 cluster, A370; it is the brightest 850 source detected in the 
survey of Small et al. (1998), and by far, the brightest SCUBA source yet detected in any of the blank field 
surveys. For all other sources only upper limits can be used to constrain the 450/xm/850/Ltm spectral index. 
Given that the noise at 450 ^m is typically 10 x that at 850 ;Um current limits on this ratio are simply not 
sensitive enough to be of much use. The only caveat might be the HDF where the exceptionally good weather 
produced a quoted limit at 450 /xm of ^ 25 mJy (Hughes et al. 1998) , implying that the two brightest 850 /xm 
sources have -§450 /-Sgso < 6 — 7; however, this sets only a weak constraint of 2; ;^ 0.3 (see Figure 6). 

K-band versus 850^m Flux 

Figure 6 shows the effects of the large negative K-correction in the submillimeter on the "submillimeter 
excess" ratio (= h'Sy{^f>^) / uSi,{2.2)). This ratio provides a useful redshift discriminant at z ^ 3, however at 
higher redshift the fact that the large positive K-corrections in both the mid-infrared and in the near-UV 
are similar means that there is little variation in the ratio over the redshift range, z ~ 3 — 6. 

Optical, Near-Infrared Colors 

One of the surprises from the early results from the Icnsed cluster survey was the finding by Small et al. 
(1998) that the optical iy — I) colors of the SCUBA counterparts (using only those sources with the highest 
probability of being correctly identified) are bluer on average than the faint field population, although the 
spread in colors is relatively large and at least a few objects appear to be extremely red. A priori one might 
have expected that the majority of the SCUBA sources, given their large infrared excess, should be extremely 
red objects (e.g. EROs: Hu and Ridgeway 1991; Cimatti et al. 1998; Dey et al. 1998). Although ULIGs are 
faint at optical wavelengths, either enough UV light escapes from the nuclear dust shroud or is still produced 
by a slowly decaying starburst in the surrounding disk. Only in the most extreme geometries is it apparently 
possible to hide all evidence of current and past starburst activity. Unfortunately, given that the optical 
colors of SCUBA sources bracket those of the faint field population, optical colors may only be of marginal 
help in identifying SCUBA counterparts. 

This ratio differs from the other template measures in that it relies on a line measurement, and thus is only 
useful once a redshift has been determined. The great majority of ULIGs (at all redshifts) appear to have 
molecular gas masses in the range -M(H2) = 0.5 — 10 x IQ^^Mq as inferred from their luminosities in one 
or more of the lower rotational transitions of CO (e.g. Sanders et al. 1991; Downes et al. 1995; Solomon 
et al. 1997). The sensitivity of current millimeterwave interferometers is sufficient to detect ^ 10^'^ Mq of 
molecular gas independent of redshift (see Radford 1994), and thus once a redshift has been determined, 
detection of the source in CO would greatly increase the probability that it indeed was the source of the 
submillimeter emission. The only drawback is that ULIGs with "warmer" than normal mid- and far-infrared 
colors can have substantially less molecular gas, e.g. M(H2) < 10^ (Evans et al. 1998), and thus fall well 
below current sensitivity limits. Therefore, given that the mid-infrared colors of SCUBA sources are not 
well-known, non-detection in CO would not necessarily rule out the identification of an otherwise reasonable 
counterpart. 

Despite the above attempts to identify the 850 //m deep-field SCUBA sources, most still do not have se- 
cure identifications. Uncertainty in the 850 /nm positions remains the major problem. Until submillimeter 
interferometers start to come on line early in the next decade, the most useful data to obtain would be inter- 
ferometer measurements at 1 mm to detect the long wavelength tail of the 850 /xm emission. Unfortunately 
the expected 1 mm fiux for most of the deep-field SCUBA sources is ^ 1 mJy, which is below the practical 
sensitivity limits of current mm-wave interferometers. Deeper "ultra-deep" VLA surveys may continue to be 
the most promising alternative. 



REDSHIFT DISTRIBUTION 



Clearly the most important unknown concerning the SCUBA sources is their redshifts. The relatively small 
fraction (~ 20-25 %) of SCUBA sources that have secure identifications has resulted in much speculation 
about the true redshift distribution. Although it is not possible to rule out distributions that peak at z ^3 
as originally advocated by Hughes et al. (1998) for the HDF, or that peak at z 1 as advocated by Bales 
et al. (1999) and Lilly et al. (1998) for the CFRS fields, there is increasing evidence that the comoving 
luminosity function of SCUBA sources peaks somewhere in between, i.e. at z ^ 1 — 3. 

The most comprehensive redshift survey of SCUBA sources is that for the 16 sources in the lensed cluster 
fields of Small et al. (1998), as recently reported in Barger et al. (1999). Of the 16 sources, two are clearly 
cluster contaminants (i.e. cD galaxies). Of the remaining 14 those with the strongest identifications are an 
interacting pair at z = 2.55 and the previously identified (Ivison et al. 1998) type-II AGN/starburst source 
at z = 2.8. The next two most reliable detections, based partly on template SEDs, are at z = 1.06 and 
z = 1.16. Eight sources are still highly uncertain, however nearly all plausible nearby counterparts have 
redshifts in the range z = 0.18 — 2.11, with two galaxies at the low end of this range likely being cluster 
contaminants. Finally, only two sources have no obvious optical counterparts, and therefore might be at 
higher redshifts (i.e. z > S). 

THE "STAR FORMATION HISTORY" OF THE UNIVERSE 

The past few years has seen a revolution in optical studies of field galaxies, which with the help of 8-10 m- 
class optical telescopes, has enabled the determination of the UV luminosity density of the Universe out to 
redshifts z ^ 4. Lookback Time (Gyr) 
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Figure 7. The global "star formation history" of the Universe {Hq = 50kms~"^Mpc~^, Qq = 0.5 is 
used for consistency with Madau et al. 1996). In the optical/near- UV, the mean co-moving rate of 
star formation is determined from the total measured rest-frame UV luminosity density of galaxies 
(solid hexagon: Trayer et al. 1998; solid circles: Lilly et al. 1996; solid squares: Connolly et al. 
1997; solid triangles: Madau et al. 1996; solid stars: Steidel et al. 1998). The optical data points 
have been corrected upward for extinction; the open symbols shown here are the correction (x4.7) 
adopted by Steidel et al. (1998). The light shaded region represents a prediction from chemical 
evolution models (Pei and Fall 1995). The dark shaded region and thick solid line represent the 
maximum contribution from far-IR/submm sources (i.e. assuming all of the far-IR/submm emission 
is powered by young stars) using models with a range of z distributions which are consistent with 
the current observations of 850 /xm SCUBA sources (Blain et al. 1998). 

Figure 7 gives the optical/near- UV determination of the mean co-moving rate of star formation over the 
redshift range z = 1 — 4 following an early analysis of galaxies in the HDF by Madau et al. (1996), who 
also noted that there seemed to be little evidence for substantial amounts of dust obscuration in the optical 
sources. It is remarkable how quickly this picture has changed. 



It is immediately obvious from the large space density of the 850 ^m sources, and their large inferred bolomet- 
ric luminosity in the far-IR/submillimcter (i.e. corresponding to ~ 30 — SOOLJ^^j), that the SCUBA sources 
could represent a significant and perhaps dominant component of the luminosity emitted by galaxies over 
cosmic time. It is also clear that a substantial fraction of the "activity" in galaxies at high redshifts (z > 1) 
is obscured by dust, and, therefore has been missed in optical/UV surveys. This is illustrated in Figure 7, 
which assuming that all of the far-IR/submm luminosity is due to star formation, plots a range of plausible 
models of SFR vs. z (Blain et al. 1998) consistent with reasonable assumptions for the redshift distribution 
of the SCUBA sources. Figure 7 suggests that the SCUBA sources dominate the optical by factors of 5-10 
at nearly all redshifts, and that their distribution with redshift is likely to be substantially different than the 
bias toward low redshifts inferred from the early optical studies. 

However, it is naive to assume that the role of dust is negligible in the optical sources, and indeed from 
studies of the optical/near- UV colors it was quickly realized that dust obscuration could have a substantial 
effect. Reddening models developed for local starburst galaxies (e.g. Meurcr et al. 1997, 1998; Calzctti 1997, 
1998) have now been used to correct the optical sources for reddening, resulting in mean model-dependent 
increases of x3-7 in the SFR (e.g. Dickinson 1998; Pettini et al. 1998; Steidel et al. 1998). Better statistical 
samples at z > 3 have also erased evidence of a turnover in the distribution at z ~ 1 (Steidel et al. 1998), 
and with an additional assumption that mean reddening corrections are relatively constant (at least over the 
observable redshift range out to z ~ 4: e.g. Steidel et al. 1998), the redshift distribution of star formation 
from the optical samples begins to resemble more closely that determined for the SCUBA sources, although 
the latter still exceed the former by factors of ~2-4 at z > 1. 

THE STARBURST-AGN CONNECTION 

Are the SCUBA sources simply powered by starbursts as is the case for the optical-UV selected samples, or 
could AGN play a significant role in producing their enormous luminosities ? 

Although relatively little is currently known about the detailed properties of the SCUBA sources, it is 
interesting to note that the bulk properties of those few sources whose positions and redshifts have been 
reliably identified with SCUBA detections are in fact quite similar to the mixture of starburst and AGN 
properties, and merger morphologies observed for ULIGs at low redshift. 

The two best studied sources from the sample of Small et al. (1997) are illustrative. SMM J02399— 0136 
at z ^ 2.8, with Ljr ^ IO^^Lq, is morphologically compact with an optical classification as a narrow- 
line "type-2" AGN (Ivison et al. 1998), and contains ^lO^^-^Mo of molecular gas (Frayer et al. 1998). 
SMM J14011+0252 at z 2.6 with Lir ~ lO^^-'^L© (Bargcr et al. 1999), is a strongly interacting/merger 
pair, with an H II- like optical spectrum and ~ lO^'^'^M© of molecular gas (Frayer et al. 1999). These two 
sources fit remarkably well into the pattern exhibited by ULIGs in the local Universe. In particular their 
molecular gas masses, optical luminosities, and optical morphologies are very similar to what is observed for 
local ULIGs (see Sanders and Mirabel 1996 for a review of local ULIGs). 

Although there are only two objects, it is interesting to note that their optical spectral types fit the trend 
with Lir that is observed for local ULIGs, as illustrated in Figure 8 which shows that the importance of 
AGN appears to increase with increasing Ljj-. At Lir > 1O^^ °~^^-^L0 10-20% of local ULIGs appear to be 
predominantly powered by an AGN with the remainder apparently powered by a starburst (although it is still 
not possible to rule out that a substantial fraction of the luminosity is produced by a heavily dust enshrouded 
AGN). At Lir > 10^"^-^ Lq nearly half of all local ULIGs appear to be powered predominantly by an AGN, 
while at the highest observed infrared luminosities (i.e. Lir ^ 10^^ L©) all of the currently identified objects 
are classified as AGN, typically as Sy2s in direct optical emission (e.g. IRAS F09105-f-4108: Kleinmann and 
Keel 1987; IRASF15307+3252: Cutri et al. 1994; IRAS F10214+4724: Rowan-Robinson et al. 1991), but 
have been shown to contain hidden broad- line regions in polarized optical light (e.g. Hines et al. 1995), or 
in direct near-infrared emission (e.g. Veilleux et al. 1997, 1998). 

Although Arp 220 (optically classified as a LINER) is often used as the template for ULIGs (see Figure 2) 
primarily due to the fact that it is the nearest such object by nearly a factor of 2.5, it is not necessarily 
typical of most UlIGs, in particular those with the largest infrared luminosities. Mrk 231 (optically classified 
as a Syl/BALQSO) may better illustrate the mixture of starburst and AGN found in many ULIGs. 
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Figure 8. The optical spectral classification of infrared galaxies versus infrared luminosity 
(Veilleux et al. 1999). 



Figure 9 shows the merger morphology apparent in deep optical images of Mrk 231, and the detection by HST 
of a large population of star clusters in the inner 1-3 kpc, presumably the remnants of a more widespread 
circumnuclear starburst whose total bolometric luminosity is currently only a very small fraction of the total 
bolometric luminosity of the system (Surace et al. 1998), nearly all of which originates from the inner few 
hundred parsecs (e.g. Matthews et al. 1987). The dense knot of molecular gas in the inner ~ Ikpc which is 
centered on the Syl nucleus, presumably represents the large reservoir of gas still available to fuel the last 
phases of the circumnuclear starburst as well as to continue to build and fuel the AGN. 




Figure 9. The advanced merger/ULIG/QSO Mrk231 - Left panel: optical image (Sanders et al. 
1987) and CO contour (Scoville et al. 1989). Right panel: i^ST B-band image and identified stellar 
clusters ('+') from Surace et al. (1998). The high resolution CO contours are from Bryant and 
Scoville (1997). 



FINAL COMMENTS 



One of the major questions at this Workshop has been the relationship of the SCUBA sources to the optically 
selected high-z population of starburst galaxies. One view is that the SCUBA sources are indeed just the 
most heavily reddened objects already contained in the optical samples. Favoring this view is the evidence 
(summarized by Steidel et al. 1998) that on average the more luminous objects in optical samples are also 
redder, such that after correction for extinction they would have intrinsic luminosities equivalent to that 
of the SCUBA sources (i.e. ^ 10^^ Lq). detect an emerged "star-formation histories" derived for the two 
populations may have similar shape might suggest that both samples are simply different manifestations of 
the same general phenomenon. 

An alternative view is that the SCUBA sources represent an inherently distinct population, for example 
the formation of massive spheroids and massive black holes, both of which are triggered by the merger 
of two large gas-rich disk galaxies (e.g. Kormendy and Sanders 1992; Kormendy and Richstone 1995; see 
also the review in Sanders and Mirabel 1996). Favoring this view is the fact that the peak in the SCUBA 
redshift distribution at z = 1 — 3 is similar to what is observed for QSOs (Schmidt et al. 1995; Warren 
et al. 1994; Shaver et al. 1998) and radio galaxies (Dunlop 1997). The relatively flat redshift distribution 
now infered for the UV starburst population at z > 1 (Steidel et al. 1998) might then better represent the 
building of gas-rich disks over a wider range of cosmic time. Also, there is no current evidence to show that 
the SCUBA detections are indeed related to the most heavily reddened optical sources, plus it is not clear 
that the far-UV luminosity of SCUBA sources is sufficient to be included in the deepest HST samples, i.e. 
L,y{2200A)/L^{80^J.m) ~ W'^-^, for Arp220-like objects (see Figure 2). 

Finally, in response to questions raised during the AGN Workshop at this COSPAR meeting, we consider 
the intruiging possibility that the large population of newly discovered SCUBA sources might be responsible 
for producing a large fraction of the X-Ray background. It is now aparent from the shape of the X-Ray 
background spectrum that most of the objects contributing to the X-Ray background are heavily absorbed 
(e.g. Fabian and Barcons 1992; Boyle et al. 1995; Almaini et al. 1998), and that they have largely been missed 
in optical surveys. Candidates for these obscured sources are the class of "narrow-line X-Ray galaxies" 
(NLXGs: Hassinger 1996), which have recently been characterized by Maiolino et al. (1998) has having 
extremely heavy obscuration along the line of sight (e.g. A^h > 10^^~^^cm~^). SCUBA sources (like local 
ULIGs) are obvious candidates for having such large absorbing columns, and if a significant fraction of the 
SCUBA sources prove to have buried AGN, then it seems reasonable to assume that they might indeed be 
responsible for producing a substantial fraction of the X-Ray background. Deep surveys with soon-to-be- 
launched X-Ray satellites (e.g. XMM, Chandra, Astro-E) could provide the answer. 

Acknowledgments. I am grateful to Karen Teramura for assistance in preparing the figures, and to 
JPL contract no. 961566 for partial financial support. 
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